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Abstract. The calculation of Euclidean distance between points is generalized to one-dimensional objects 
such as strings or polymers. Necessary and sufficient conditions for the minimal transformation between 
two polymer configurations are derived. Transformations consist of piecewise rotations and translations 
subject to Weierstrass-Erdmann corner conditions. Numerous examples are given for the special cases of 
one and two links. The transition to a large number of links is investigated, where the distance converges 
to the polymer length times the mean root square distance (MRSD) between polymer configurations, 
assuming curvature and non-crossing constraints can be neglected. Applications of this metric to protein 
folding are investigated. Potential applications are also discussed for structural alignment problems such as 
pharmacophore identification, and inverse kinematic problems in motor learning and control. 
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1. Introduction 

The standard variational definition of distance can be generalized to higher dimensional objects such 
as strings or membranes. In a previous paper [1], one of us has introduced the formalism for this 
calculation. Consider first zero-dimensional objects (points). The distance between two points A and 
B is defined through a transformation that takes A to B, an object of dimension one higher than the points 
themselves (here one-dimensional). The transformation minimizing the arc-length travelled between A 
and B gives the scalar distance V*. The differential increment of arc-length may be defined as either 
a/1 + (dy/dx) 2 + (dz/dx) 2 dx, or without the assumption that y,z are functions of x, parametrically. 
To be specific, introduce a "time" parameter t such that < t < T, and r(0) = r A , r(T) = r B , and 
r(f) = (x(t), y(t),z(t)). The distance between r A and r B can be found variationally 0: 

V* = V [r* (t)] where r* (t) satisfies (la) 

5 fdt{g^{t)x u {t)) l/2 = 0. (lb) 
Jo 

or 5 dtVt 2 = (Euclidean metric) (lc) 



o 



Here we have let x = dx/dt, and r = dr/dt. The boundary conditions on the extremal path are r* (0) = r A 
andr*(r) =r B . 

Taking the functional derivative in eq. (PTcl) gives Euler-Lagrange (EL) equations for the Lagrangian 
£ = Vi 2 ~: 

dt\dr ) 

or v =0 (2) 

with v the unit vector in the direction of the velocity. 

Since the derivative of a unit vector is always orthogonal to that vector, equation © says that the 
direction of the velocity cannot change, and therefore straight line motion results. Applying the boundary 
conditions gives v = (r B — r A )/ |r B — r A |. However, any function \{t) = \v a (t)\ \ satisfying the boundary 

conditions is a solution, so long as J Q dt \v (t) \ = |r B — r A |. The solution is reparameterization-invariant. 
Then the extremal functional r* (t) is given by 



r*(0=r A + -^-^ / dt\ Vo (t)\ (3) 
l r B — r A | Jo 

and the distance by 

V*= [ dtV^= [ dt |v„(0 1 = |r B — r A | (4) 
Jo Jo 

which represents the diagonal of a hypercube, as expected. At this point we could fix the parameterization 
by choosing \v (t) \ = |r B — r A | /T (constant speed), for example. 

The extremal transformation © is also a minimum. In section 12.41 we will give the sufficient 
conditions for an extremum to be a (local) minimum, where we will return to this example. 

The above idea can be generalized to space curves, surfaces, or higher dimensional manifolds [QQ|. 
The distance is defined through the transformation between the objects that minimizes the cumulative 
amount of arc-length travelled by all parts of the manifold. 



2. Distance for polymers or strings 



Describing the transformation r(sj) between two space curves r A (^) and r B (s) requires two scalar 
parameters: s the arc-length along the space curve, and t the "time" as in the above zero-dimensional 
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case measuring progress during the transformation. The boundary conditions are then r(s, 0) = r A (s) and 
r(s, T) = r B (s). The minimal transformation r*(s,t) is an object of dimension one higher than A or B, i.e. 
it yields a distance that is two-dimensional. 

The distance V* =V[r*(s 1 t)], where the functional V[r] is given by 



V[r] = ds dt Vr 2 . (5) 



Here we have used the shorthand r = r(s,t) = (x(s 1 t) 1 y(s 1 t),z(s 1 t)) (a 3-vector), and r = dr/dt. 

It has been shown previously that the problem of distance does not map to a simple soap film, nor to 
the minimal area of a world-sheet (which corresponds to the action of a classical relativistic string) [0Q. 

Formulated as above, the string can contract and expand arbitrarily in order to minimize the distance 
travelled. The transforming object is akin to a rubber band, and all points on r A (s) will move in straight 
lines to their partner points on r B (s) to minimize the distance. It is worth mentioning that protein chains 
for example only change their length by about one percent at biological temperatures. To accurately 
represent the transformation of a non-extensible string, a Lagrange multiplier X(s,t) must be introduced 
into the effective Lagrangian, weighting the constraint: 



r' 2 = 1 , (6) 
where r' = dr/ds. 

Under this constraint, points along the string can no longer move independently of each other, but 
must always be a fixed (infinitesimal) distance apart. The tangent vector t = r' is now a unit vector, and 
the total length of the string is L = J^ds x/r' 2 = J^ds. 

Consider the minimal distance transformation between two configurations r A (s) and r B (^) of an ideal 
polymer of length L. Let us derive the EL equations for this case. From equations © and ©, the 
effective action is 



V =11 dsdtC(t,r') 

Jo Jo 



(7a) 

where C = Vr 2 - 1 ( Vr' 2 - 1 ) (7b) 



and the Lagrange multiplier 1 = l(s, t) is a function of both s and t. The extrema of the distance functional 
V in (TTab are found from 5V = 0. Taking the functional derivative gives EL equations 0]|: 

v = lK + l't. (8) 

where v is the unit velocity vector, t is the unit tangent vector, and n is the curvature vector. In eq. © 
we see explicitly that if the non-extensibility constraint is set to zero (1 = 0), all points on r A (s) move in 
straight lines to r B (^). 



2.1. Discrete chains 

To make the problem more amenable to solution, we can discretize the spatial variables while letting 
the time variable remain continuous, i.e. we implement the method of lines to solve eq. ©. Rather 
than directly discretizing eq. ® however, it is more natural to consider a discretized chain as shown in 
figure Q] from the outset, and to calculate the EL equations for this system. This recipe then gives the 
same result as properly discretizing eq. ©. For the discretized chain, the constraint in eq. © becomes 
| Ar| = As = L/(N — 1), giving the length of each link. As the number of beads N — > oo the system 
approaches a continuous chain. For finite Af, the Lagrangian becomes a function of the positions and 
velocities {r,,r,} of all beads i, 1 < i < N+ 1. We use the shorthand notation £(r,-, r,-). 

This recipe yields the distance metric for an ideal, freely-jointed chain, which has no non-local 
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Figure 1. Continuum (a) and discretized (b) polymer chain. The EL equation for the continuum polymer is 
a nonlinear (vector) PDE, while the EL equations for the discretized polymer are a set of nonlinear ODEs. 

interactions and no curvature constraints. While this approximation is often used as a first step, real 
chains may behave quite differently for several reasons. In many cases, the configuration which is an 
energetic minimum is a straight line, or a single conformation dictated by the chemistry of the polymeric 
bonds. At finite temperature, energy in the bath induces conformational fluctuations. Real polymers 
also cannot cross themselves, and because of their stereochemistry also take up volume. We leave these 
interesting features for later analysis. 

Equation © for the discretized chain becomes N constraint equations added to the effective 
Lagrangian: 

N 

J^Xv+i (V(r ; - + i-r ; -) 2 -A/) 
i=i 

where each A^+i = A ; ,r+i (0 is a function of t, and t^w+i — 0. Letting 1 = 21 As and r,- +1 /,- = r, +1 — r,- we 
rewrite this strictly for convenience as 

Eh,i+l ( r i+l/i 
~T~ I As 2 

We next convert to dimensionless variables by letting r = ( As)r. To simplify the notation, from here 
on we simply refer to r as r. The distance for the discretized chain becomes 

V[r h r i ] = As 2 [ dt (9) 
Jo 

with effective Lagrangian 

The derivatives r and r, +1 / ( are raised to different powers in ( flOl ), however so long as r,- +1 /,- satisfies 
the constraint |r(+i/;| = 1, the EL equations for r f (f) will be the same whether the constraint ^r 2 +l ^ = 1 
or rf+i/i = 1 is used. 

The reparameterization invariance present for point particles (c.f. section [B is still present for 
beads on the chain, but the parameterization of arclength along the chain is taken to be fixed by the 
discretization. 



N / / f ^ 

(10) 
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Figure 2. General variations of a functional with fixed end points allow for broken extremals. In the text 
we derive the extra "corner" conditions for a piecewise continuous path to still be extremal for our distance 
functional. 

2.2. General variation of the distance functional 

For reasons that will become clear as we progress, we consider the general variation of the functional 
V, allowing for broken extremals. That is, we allow the curves describing the particle trajectories to be 
non-smooth in principle at one or more points in time. Consider the case of one such point at time t\ . 
The distance can be written as 

V= [dt £(r h r i ) + [ dtC(Ti,Tt) (11) 



The space-trajectories of the particles must be continuous at time t u so rfa — e) = + e), or in 
shorthand: 

r / (fr)=r i (f+) . (12) 

Let r,-(f) and ?,•(*) be two neighboring trajectories from r,(0) = r A( - to r/(r) = r B; (see figure |2]). 
Neighboring curves will differ by the first order quantity h,-(f) = r,-(f) — r,-(f). The fixed boundary 
conditions at t = 0, T dictate that h ( (0) = h,-(T) = 0. The difference in distance between the two 
trajectories is 

AX) = X)[r i + h < ]-X)[r i ] 

rh+Stj rt\ rT i-T 

= / dt C{Ti+hi,Ti+hi) - / dtC(r U Ti)+ / J^^ + h^r^ + h,) - / dt£(rQQ) 

JO J() Jti+Sti Jti 

Taylor expanding the Lagrangian to first order in h,:* 

N 

£^£(r,-,r,) + ^(£ r ,-h, + £ r ,-h,) 

i=\ 

and integrating by parts using the fixed boundary conditions at t = 0, T, the difference in distance up to 
first order in h; is 

AD « \\ ^ (c ri - | A ) ■ h, + [d, j: - • h, 

+ £((,-)&,-£(<+)&,+ A-H,-- ^A-Hf < 14 > 

i i 

| i We use the notation F r = dF/ dr, Fi = dF/ dr. 
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with the shorthand C{t) = £(r,-(f),r,-(f)). 



2.3. Conditions for an extremum 



The variation 5V differs from AD above only by second order terms. Then for the transformation from 
{r A , } to {r B , } to be an extremum, 8V = 0. Thus, the EL equations (in the top line of eq. ([14)) ') must vanish 
in each regime [0, ti), (h,T]. Using the form of the Lagrangian in eq. (flOl) , the EL equations become: 



Vi+li 2 r 2 /i =0 

V 2 -ll2 1*2/1 +1231*3/2 = 



(15a) 
(15b) 



Vw-liv-i,jvrjv/(iv-i) = (15c) 

According to equation (fl4l) there are additional conditions for the transformation to be an extremum. 
To find these first note that up to first order (see figure O 

h i (t 1 )^5r i ( y t l )-r i (t l )5t l . (16) 

Then the first variation in the distance is 



SV 



+ 




8t\ 



<*r,-(/i) 



(17) 



which must vanish at an extremum. Because the variations Sr { and 8t\ are all independent, the terms 
in square brackets in equation (TT71) must vanish. Writing these expressions in terms of the conjugate 
momenta p, = and Hamiltonian Ti = 2~2i*i • P; — £ gives the conditions: 

(18a) 



P/l 

i 

Til 



P/l 

t 

n\ 



(18b) 



These conditions are called the Weierstrass-Erdmann conditions or corner conditions in the calculus of 
variations 0. 

According to the Lagrangian in equation (flOl) , the Hamiltonian is given by 



U 



-E^( r W-l) 



;=1 



which is identically zero, so corner condition (|18b| ) provides no further information. 
The conjugate momenta according to (flOl) are given by 



(19) 



Therefore, according to corner condition (I18al) . extremal trajectories cannot suddenly change direction: 
each Ti(t) follows a smooth path continuous up to first derivatives in the spatial coordinates. 

The fact that one corner condition provided no information due to the vanishing of the Hamiltonian 
is related to our choice of parameterization in formulating the problem. For example, in the case of the 
distance of the single point particle mentioned in the introduction, the Lagrangian may be defined either 

through independent variable x as £W = a/ 1 + y' 2 + z' 2 (with e.g. y' = dy/dx), or parametrically through 

independent variable t as = yf&. The conjugate momenta are then either =y' 



i/2 . 
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and = z'/^l+y^ + z' 2 , or Cf = r/|r| = v. The Hamiltonia are either H {x) = \/^\+y' 2 + z' 2 
or Ti.^ = — r • (r/|r|) = 0. The corner conditions can be shown to be equivalent for both choices 
of independent variable: for £W they give v(?f ) = v(/j + ), so that the direction of the tangent to the 
curve cannot have a discontinuity. Together, the Hamiltonian and two conjugate momenta for £w can be 
interpreted as components of the unit tangent vector to the curve, i.e. t(x) = (i+j'j+z'k) / a/1 +y' 2 + z' 2 , 
and so once again the corner conditions enforce a continuous tangent vector, here t(xf) = t(x+). 

2.3.1. Boundary conditions In the continuum limit, the boundary conditions on r(s,t) are r(s,0) = 
r A (s), r(s,T) = r B (s) where r A and r B are the two configurations of the polymer. For discrete chains, 
these boundary conditions become 

{r,(0)} ={rf} (20a) 

{r,(r)}={r| B) }. (20b) 

There are also boundary conditions that hold for the end points of the chain at all times. From 
equations (I15a[ I15cl) we see that there are three solutions for the end points of the chain: 

1) If 1 7^ 0, purely rotational motion results. This can be seen by taking the dot product of eq. (I15al) 
with Vi, which yields l 12 Vi ■ r 2 /i = 0, so the velocity of the end point is orthogonal to the link. The 
rotation must be about a point that is internal to the link, i.e. on the line between points 1 and 2 for end 
point 1 . This can be seen straightforwardly for the case of one link by removing point 3 from equations 
(1 1 Sab and (|15b| ). Then the accelerations v, must be in opposite directions. This can only occur if rotation 
is about a point on the line between points 1 and 2. 

2) If 1 = 0, Vj = 0, and straight-line motion of the end point results. 

3) Writing out the time-derivative in (1 1 5ab yields 

v?vi - (vi -vi)vi = — 112 | Vi | 3 r 2 /i (21) 
which has the trivial solution Vi = 0. The end point can be at rest, while other parts of the chain move. 



2.4. Sufficient conditions for a minimum 

For a transformation to be minimal, it is necessary, but not sufficient, that it be an extremum. We now 
derive the sufficient conditions for a given transformation to minimize the functional ©. We describe 
the formalism in some detail because it is not typically taught to physicists- for further reading see for 
example reference [2J. This section can be read independently of the others, and might be skipped on 
first reading. 

According to Sylvester's criterion, a quadratic form J2ij^-ij x i x j 1S positive definite if and only if all 
descending principle minors of the matrix ||Ay|| are positive, i.e. 

>0, ... ,det||A y || >0, (22) 

and a function F of x = (xi,X2, ■ ■ ■ ,x n ) has a minimum at x* if the Jacobian matrix \\d 2 F /dxidxj\\ is 
positive definite at the position of the extremum (where OF /dxj = 0). 

For a function to be a minimum of a given functional, it must satisfy similar sufficient conditions. 
Consider again the difference in distance between two trajectories in Taylor expanding the 

| * We ignore corner conditions for purposes of the derivation. It can be shown that they do not modify the result. 



A,, >0. 



An A i2 
A21 A22 



>0, 



An 


A 12 


A13 


A21 


A 22 


A 23 


A31 


A 32 


A33 
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Lagrangian to second order in h, : 

AV = V[r i + h i ]-V[r i ] 

r 

dt£(Ti+bi,Ti+hi) 



dtC(r h Ti) 



r T N i 3N 

dt ^2 (A; ■ h,- + C fi ■ hi) + - ^2 {ZxiXjhihj + 2C X , X //,/// + £ V v,////i/ 



(23) 



At an extremum, the first order term in (|23T) is zero, and A£> « 5 2 £>, the second variation. For the 
extremum to be a minimum, 5 2 V > 0. From eq. (flOl) . the matrix ||A,.vJ = ||0||. Assuming ||£r,-.v,|| is 
in general a symmetric matrix, i.e. C XiX - = C X Xi , the second term in the quadratic form of (|23T) may be 
integrated by parts to give: 

1 f T r • • 1 

' (h\m) + (h\Qh) , (24) 



5 2 P 



c/z 1 



where we have let \h) denote the vector (h\,fi2, ■ ■ ■ ,hw), and used the shorthand P and Q for the matrices: 



\\£- ■ II 

*^(\\c \\--\\c-\\\ 

From (fTOl) the explicit form for these matrices may be calculated. P is block diagonal: 



Q(0 = HQ/. 



(25) 



'j 




• 


■ 


/( 2) • 

<J 


• 


■ 


ij - 



(26) 



with each block matrix having elements 

\\tM\\ = J_^..r^ 2 -jc^i c 
II 1 '; II I -/i3 ^°'v r x i x 7 



y 2 + z' 2 -xy' 
— xy x 2 + i : 
— xz 



—xz 
-yz 
yz x 2 +y 2 



(27) 



(particle J) 



Interestingly the numerator of (1271 has the form of an inertia tensor for a point particle in velocity-space. 
The matrix Q is block tri-diagonal, because the spatial derivatives in (|25T) couple each bead to its two 
neighbors. Using indices I, J to enumerate beads and i, j to enumerate x, y, z components for each bead: 

\\Qjj,ij\\ = % {Sij — Sij-i) +h,j+i (5ij — Sij+i)] 
"li 2 l -I12I 

-ll2l (ll2+l 23 )l -I23I 

-1 23 1 (1 23 + 134) 1 -I34I 



or 



Q 



-l 



Ar-lJV 



1 1 



Af-lJV-l 



(28) 



For the transformation r*(?) to be a minimum of V[r], the functional (1241) must be positive definite 
for all To derive the conditions for this, we can temporarily ignore the fact that (1241) arose from the 
second variation of ©, and treat (|24l) as a new functional of the function \h(t)) = \hi(t), . . .h 3N (t)). We 
then ask what \h(t)) extremizes (|24|) . If 5 2 V > we expect that the only extremal solution would be the 
trivial one: \h(t)) = |0), at least for small variations of the hi(t). That is, changing the transformation 
{r*(t)} from that which extremized © to a neighboring transformation {r*(t) +h,(?)} would increase 
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the distance travelled. 

The system of 3iV EL equations for \h) from ((24)) is 

-jL|pji) + |Qft) = |0) (29) 

with boundary conditions 

\h{0)) = \h{T)) = |0) . (30) 

Equation (|29l) is referred to as the Jacobi equation in the calculus of variations. 

First note that if \h) satisfies the system of equations in (1291) as well as the boundary conditions (1301) . 
then integration by parts gives 

S 2 V = J^dt ((h\Ph) + (h\Qh)^ = ^dt(h\-j t (ph^+Qh)=0. (31) 

This means that for 5 2 V to be > 0, any nontrivial \h(t)) which satisfies the boundary conditions must 
not itself be an extremal solution of the Jacobi equation, otherwise solutions |r* (?) ) perturbed by any 
constant times | h(t)) are themselves extremals. One may think of this by analogy as the necessity for 
the absence of any "Goldstone modes", where excitations by various C\h{t)) would lead to a family of 
curves with zero cost in action, and thus zero effective restoring force, between them. 

Alternatively we can ask what equation h = \h) must satisfy if the EL equations are satisfied for both 
£(r,r) and the neighboring extremal C(r + h,r + h). Taylor expanding C(r + h,r + h) in 

£ r (r + h,r + h)-4-£r(r + h,r + h) =o 

dt 



gives 



-!(/*-h) + (£.-^)-k = 



which is exactly Jacobi's equation (1291) with definitions (1251) . 

From here on, it is much simpler to elucidate the central concepts for sufficient conditions using the 
case of a single scalar function h(t). The analysis can be generalized to the multi-dimensional case with 
a bit more effort, but the conclusions are essentially the same and so they will simply be stated along 
with the conclusions for the ' 1-D' case. For further details see 0. 

We write equation (l24l) in 1-D as: 



T 

dt(Ph 2 + Qh 2 ) (32) 



It was realized originally by Legendre that the integral could be brought to simpler form by adding zero 
to it in the form of a total derivative. Since 



/ dt^- (w(t)h 2 ) = 
Jo dt 



for any w(t) so long as h{t) satisfies the boundary conditions (l30l) . we can add it to the integral in ([32 
and seek a function w(t) such that the expression 



5 2 V = i f dt (Ph 2 + 2whh + (Q + w) ti 
2 Jo ^ 
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may be written as a perfect square. This yields the differential equation 

P(Q + w) = w 2 (33) 

for w(t), and second variation 

S 2 V[h] = i [ dtp(h + -h) 2 . (34) 



2J \ P 

Therefore a necessary condition for a minimum is for P > 0. The analogous condition in the multi- 
dimensional case is for the matrix ||P|| to be positive definite. 

If the differential term h + jh in (|34l) were equal to zero for some h(t), the boundary condition 

h(0) = would then imply h(0) = and thus h(t) = for all t by the uniqueness theorem as applied to 
this first order differential equation. 

Therefore the functional (l34l) is positive definite if, and only if, 

1. ) P > , 

2. ) A solution for eq. (1331 exists for the whole interval [0, T]. 

In general, there is no guarantee of condition (2) even if condition (1) is valid. For example if P = 1 , 
Q = — 1, (|33l) has solution w(t) — tan(? + c), which has no finite solution if \T\ > n. t 

If (|33l has a pole at say t, then for the integral (13~4l ) to remain finite, — > 0. This point is said to be 
conjugate to the point t = 0, i.e. it is a conjugate point. 

Moreover, equation (1331 is a Riccati equation, which may be brought to linear form by the 
transformation w(t) = —PH/H, with H{t) an unknown function. Substitution in (1331 gives 

-j t (PH)+QH = (35) 

which is precisely equation (|29l )- the Jacobi equation for h(t). 

This means that for equation (1331) to have a solution on [0, T], H(t), as given by the solution to (T35l) . 
must have no roots on [0, T]. But because equation (1331) holds for h(t) as well, h{t) must have no roots 
(conjugate points) on [0,T]. Because /z(0) = h{T) = 0, the only way to extremize d32l) is to satisfy 
eq. (1331 ) with the trivial solution h(t) = 0. If ^ for < t < T then it would mean that there was a 
conjugate point at t = T. 

In the multi-dimensional case an extremal \h) is one of 3./V vectors satisfying equations (|29l) , i.e. 

\h^) = \h^ . . .hyj), 1 < a < 3N. A conjugate point is defined as a point where the determinant 
vanishes 



det 



o 



The conditions for a transformation to be minimal are then: 

1. ) The transformation \r*(t)) = {r*(t)} is extremal, 

2. ) Along \r*(t)), the matrix P(f) = L iiij is positive definite, and 

3. ) The interval [0, T] contains no conjugate points to t = 0. 
The above ideas can be made clear with a few examples below. 



2.4.1. Distance between points From the effective Lagrangian £ = Vr^, P = is given in 

equation (1271) , which has determinant detP = 0, and so is not positive definite. This is due to our choice 
of parameterization. If we break symmetry by choosing one spatial direction as the independent variable, 

| t Because reparameterization invariance in our problem, the value of T is adjustable, however precisely because of this 
invariance, det ||P| = and so is no longer positive definite. We discuss this problem and its resolution below. 
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C (x,y',z') = a/ 1 + y' 2 + z' 2 (with e.g. / = dy/dx andx < x < xi). Then 

P _ 1 fl+z' 2 -y'z'\ 

with positive definite determinant det||P|| = (l +y' 2 + z /2 ) > for any trajectory. From eq (|25l) . 
||Q(0II — l|0||- Along the extremal, where y(x) = ax + y , z{x) = bx + zo, equation (l29l) gives P h' = c, 
with c a constant vector and P a positive definite matrix of constant values with respect to x. Solving this 
first-order equation gives straight line solutions for h(x). Because h(x ) = 0, there can be no conjugate 
points, and because h(xi ) = 0, the only solution to d29l) is the trivial one, and the extremum is a minimum. 

2.4.2. Geodesies on the surface of a sphere Taking the azimuthal angle <p as the independent variable, 
and polar angle 9(4>) as the dependent variable, the arc-length on the surface of a unit sphere may be 
written as 

V[6]= I d<\> v 7 ' 2 + sin 2 . (36) 

J <f) 

The EL equations give the extremal trajectory as cos = Asin cos0 + 5sin sin0 with A, B constants. This 
is the equation of a plane z = Ax + By, which intersects the surface of the sphere to make a great 

circle. The scalar P = C" = sin 2 / ( /2 + sin 2 ) 372 which is always positive. To simplify the problem, 
let </>o = 0, and (</> ) = (4>i) = n/2, so the great circle lies in the z = plane. Along this extremal P is 
constant and equal to 1, while Q = — 1. The second variation, eq. (1321) . is then (1/2) J^dcj) (h 12 — h 2 ). 
The corresponding Jacobi equation, h" + h = 0, must not have a root between [O,0i]. The nontrivial 
solution to the Jacobi equation satisfying the initial condition h(0) = is h(<p) = Csin0, which has a 
conjugate point at <p = 7T- Thus for the extremal curve to be minimal, (j>\ must be < n, the location of 
the opposite pole on the sphere. If 0i < n, there is no extremal solution for h((p) other than the trivial 
one which satisfies the boundary conditions. It is instructive to look at the arc-length under sinusoidal 
variations around the extremal path which satisfy the boundary conditions h(0) = /z(0i) = 0, so that 
9(4>) = ir/2 + h((fi) = 7r/2 + esin (7r0/0i) . Inserting this into eq (|36l) above and expanding to second 
order in e, we see that first order terms in e vanish, and the difference in distance from the extremal path 
is AD = (e 2 /40i) (71" 2 — 2 ). For <p\ < n this is always greater than zero indicating the extremal is a 
minimum. For (fii > n this is always less than zero indicating the extremal is a maximum with respect to 
these perturbations: the length may be shortened. When <pi = ir, AD = to second order. When h((p) 
represents the difference between great circles AV is precisely zero. 

2.4.3. Harmonic oscillator It is not widely appreciated that the classical action for a simple harmonic 
oscillator is not always a minimum, and indeed in many cases can be a maximum with respect to 
some perturbations. The action for a harmonic oscillator with given spring constant is proportional to 

S[x] = J^dt i(i 2 —x 2 ), which has EL equation x + x = 0. Taking the specific initial conditions x(0) = 1, 
x(0) = 0, the extremal solution is x(t) = cost. The scalar P(t) = Ca = 1, which is always positive 
and satisfies the necessary conditions for a minimum. The scalar Q = £ xv — jf t C X x = — 1- The second 

variation 5 2 S[h] = \ J^dt (h 2 — h 2 ), which has Jacobi equation h + h = 0. This is the same Jacobi equation 
as that for geodesies on a sphere, so the sufficient conditions will parallel those above. The boundary 

condition /i(0) = gives h(t) = A s'mt, with conjugate points at t = nir, n= 1,2, This means that the 

action is a minimum only so long as T < ir, i.e. a half-period. If we let x(t) be the extremal solution 
plus a sin perturbation satisfying the Jacobi equation at the conjugate points: x(t) = cost + es'mt, then 
the difference in action from the extremal path becomes AS = (e 2 /4T)(Tc 2 — T 2 ). This result is exact 
because the action for the oscillator is quadratic (as opposed to the action for geodesies). When T < ir, 
AS > indicating the extremal is a minimum. When T is larger than a half-period, AS < and the 
extremal trajectory is a maximum (with respect to half-wavelength sinusoidal perturbations), and when 
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a b c 



Figure 3. Possible (a,b) and impossible (c) straight line transformations between links AB and A'B 1 . Figure 
b shows a straight line transformation where the initial and final states do not lie in the same plane. In the 
text we derive the conditions for the possibility of a straight line transformation between links. 

T = it, the end point is the conjugate point and AS = 0. 

We discuss sufficient conditions further below in the context of minimal transformations for links. 



3. Single Links 

In the limit of one link, equations (|15aH15c|) reduce to: 
v A + lr B / A = 

v B -lr B/A = (37) 

where we have let A represent point 1, B point 2, and 1 = \ n . The link has length 1 in our dimensionless 
formulation, so the vector r B/A could also have been written as a unit vector r B/A . 

Both points A and B are end points and satisfy the boundary conditions of section |2~.3. II This means 
that points A and B move by either pure rotation, straight-line translation, or remain at rest. The initial 
and final conditions may be written r A (0) = A, r B (0) = B, r A (T) = A', t b (T) = B'. 

The link in our problem has direction, so A must transform to A' and B to B' . We will often use 
arrowheads in figures to denote this direction. 



3.1. Straight line transformations 

As a first example, consider the two links shown in figure [3^. The four points A, 5, A', B' need not lie in a 
plane (see for example fig[3j5). Let angle ZBAA' = a be obtuse. We draw straight lines from A to A' and 
B to B', and ask whether such a transformation is possible. We can thus derive the following rule: 
• For a straight line transformation to exist between two links, opposite angles of the quadrilateral made 
by AB, A'B', AA, B~W must be obtuse. 

Let the length that point A travels be x A , i.e. we imagine the point A' and the distance x A = \AA'\ to be 
variable. The length r B that point B travels is then a function of x A and the original angle a, r B (x Al a). We 
can now find conditions on the angle b = ZBB'A' such that the transformation is possible. 

After some distance x A travelled by point A, the length of the line from B to A' is 

BA' = x 2 A + 1 — 2x A cos a 
= r\ + 1 — 2r B cos b 

so that 

r B (x A , a) = cos b ± \/cos 2 b + f(x Al a) 
with f(x Al a) = x 2 A — 2x A cosa. Since a is obtuse, / > when x A > 0, and so the positive root must be 
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A B 

(A) An example of a set of link configurations connected by a straight-line transformation. The link rotates 
clockwise as it translates to allow the end points to move in straight lines. The translation can proceed no farther 

than the end points AB and A'B', which have link vectors AB or A'B' that are perpendicular to one or the other of 
the vectors v A or v B . The totality of states thus connected forms a "bowtie". (B) A bowtie where the terminal 

states AB and A'B' happen to cross each other. 

taken for r B to positive. When x A = 0, /(0, a) = 0, and 

r B (0,a) = cos&+ |cos£>| = 

Therefore b must also be an obtuse angle. If two opposite angles are obtuse, then the other two angles 
must be acute. This concludes the proof that the above conditions are sufficient. An additional proof that 



they are necessary is given in Appendix A 

We readily see that figure HA is one pair of a larger set of straight line transformations that can 
continue until one or both of the obtuse angles reaches 90°. This collection forms a "bow tie" of 
admissible configurations, as in figure HI Note that straight lines in the quadrilateral may cross as in 
the transformation from A, B to A', B' in figure HI Trivial translations of the link without any concurrent 
rotation are a special case of general straight line transformations. 



3.2. Piece-wise extremal transformations: transformations with rotations 

An immediate question is the nature of the transformation between AB and A'B' in figure [3]C, where 
opposite angles of the quadrilateral are not obtuse. Recall our link has direction so A cannot transform 
to B' . Then direct straight-line solution is not possible due to the constraint of constant link length. 

The only remaining solution is for the link to rotate as part of the transformation. Consider first the 
rotation of link AB. The EL equations (1371) allow for pure rotations about A, B, or a common center along 
the link. Likewise for link A'B'. 

The rotation can occur from either link AB (fig [5^) or link A'B' (fig [5p). After the link rotates to a 
critical angle, it can then travel in a straight line. The extremals are broken in that they involve matching 
up a piece consisting of pure rotation with a piece consisting of pure translation of the end points of the 
link. Where the pieces match they must satisfy the corner conditions (|18alll8bT) . This means that the end 
points cannot suddenly change direction, a situation which is only satisfied by a straight line trajectory 
that lies tangent to the circle of rotation. 

From figure HI we see that a straight line transformation exists only when an angle between a link and 
one of the straight line trajectories reaches ir/2. 

The critical angle that link AB must rotate is then determined by the point where a line drawn from B' 
is just tangent to the unit sphere centered at point A, point B{ in figure [5^. There is generally a different 
critical angle if the rotation occurs at link A'B' as in fig 0B. It is shown in | Appendix B that in general the 



critical angle is determined by drawing the tangent to a circle or sphere about one of the link ends. 

If the rotation was about a common center, we see that one or another of the link ends would violate 
a corner condition, so the rotation must be about one of the link ends. 
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Figure 5. Transformations between two links involving broken extremals consisting of rotation and 
translation, (b) is the global minimum, with shortest distance travelled during the transformation, (a), 

(c), and (d) are local minima, (e) is extremal, but not minimal as the trajectory of arc B'B\ passes through 
a conjugate point. 

According to eqs. (1261) and (1271) . the matrix P has a determinant of zero due to the parametric 
formulation in the problem and so is not positive definite. To show that the transformations in fig.[5^,b 
are indeed minimal, we need to then express the problem in non-parametric form. To do this, let the 
independent variable be the angle of the link with the vertical. Then the displacement x along the line 
AA' is the unknown function of to be determined by minimizing the total arc length travelled. This 
distance can be written as 



V[x] 



V 2 + 2x'cos + 1 



In this formulation, the scalar quantity P() = C x i x i becomes 



P() 



sin 



(x /2 + 2^cos + 1) 



3/2 



which is always > except for the isolated point = 0, in particular it is positive along the extremal 
trajectory which is necessary for a minimum. So we conclude that the transformation with the smaller 
angle of rotation in fig [5t> is here the global minimum, and the other transformation (fig [5k) is a local 
minimum. 

Figure [5^ is also an extremal trajectory, satisfying corner conditions, and with positive definite P. 
However it is not a local minimum because the trajectory passes through a conjugate point (denoted by 
point CP, where the dotted line along A'B' meets the great circle about A'). According to the results in 
section [2~.4.2l if the extremal trajectory (a great circle) traverses an angle larger than n radians, it passes 
through a conjugate point and thus becomes unstable to long-wavelength perturbations. Transformations 
involving rotations about points B or B' in figure |5]both have conjugate points and so are not minimal. 
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The transformation in fig. [5fc does not pass through a conjugate point and so is in fact another local 
minimum. The part of the extremum along the straight line section of the trajectory has no conjugate 
points as discussed above. 

3.3. Systematically exploring transformations by varying link positions 

We can investigate what happens to the minimal transformation when one of the link positions or angles 
is varied with respect to the other. Let us start by putting the two links head to tail as shown in figure [6]A 
The distance between them is 2 by simple translation of link end points. 

We can now increase the angle between the two vectors by rotating the right link for example, as in 
figures 0B-H. So long as the angle between the two vectors is less than 90°, one link may slide along 
another and the distance is unchanged (figs[6j\-C). This is a special case of the transformations shown in 
figure H] (compare for example figure [6j3 with the middle three unlabelled links in that figure). 

Beyond 90° however, the transformation must include rotation. Fig [60 has an angle of 150°. The 
minimal transformation first rotates, for example with the tail of the horizontal black arrow fixed, and 
the head tracing out the blue arc, until the critical angle is reached, where a straight line made from the 
final arrowhead (at the top of the figure) is just tangent to the circle made by the blue arc. This state 
is indicated by a red link in figure [6j3. The link then translates to its reciprocal position at the opposite 
end of the bowtie, denoted by a second red link (c.f. also figure @JB). At this point the arrowhead has 
completed the transformation. Finally the tail rotates into its final position. The total distance travelled 
is slightly larger than 2. 

When the angle between the vectors is 120° as shown in [6ji, the transformation consists of pure 
rotations. Taking the initial state to be the horizontal black vector, the link first rotates about its fixed tail, 
the head tracing out the blue arc, until the link reaches the state shown in red, where the position of the 
arrowhead has reached its final end point. Then the link rotates about its head until the position of the 
tail reaches the final state. 

When the angle between the links is larger than 120° as shown in figs[6f'-G, the transformation must 
involve rotation about an internal point along the link. Let points A and B denote the tail and head of 
the link respectively. If an infinitesimal rotation A occurs about an internal point P, the increment in 
distance travelled is 

AV= |r B/P |A + |r B/A |A = A 

which is independent of the position of the instantaneous center of rotation (ICR). This means that there 
are an infinity of transformations all giving the same distance, depending on the time-dependence of 
the ICR. Two simple alternatives with only two discrete positions of ICR are shown in figures [6f,G. 
Specifically, in figure [6f, the horizontal black vector first rotates about its tail to the red configuration, 
which is a mirror image of the final black vector. Then rotation is about an internal point determined by 
the intercept of the red vector with the final black vector, with end points tracing out the green arcs. In 
figure [6j3 the two ICRs are both internal and determined by the intercepts of the initial and final states 
with the red vector shown. 

Figure [6)3 depicts the transformation for overlapping, opposite pointing vectors. Rotation can now 
only occur about one point in the center of the vectors. 

Figure |7] illustrates what happens when one of the links is translated with respect to another, starting 
from two different scenarios shown in |7J\ and UB- In EK, the tail of the vertical link is displaced 
(1/3,-1/3) with respect to the tail of the horizontal link. The minimal transformation is a pure rotation 
by 7r/2. 

In figure [TJB, the tail of the vertical link is now displaced to (2/3, —1/3). Pure rotations again give a 
distance of w/2. Rotation about a point on the horizontal link that is equidistant from both arrowheads 
transforms the initial arrowhead to the final (red intermediate state). Then rotation of the tail about the 
arrowhead transforms to the final state. 

In figure the minimal transformation first involves a translation by sliding the arrowhead along 
the vertical, until the arrowheads overlap (red intermediate state). The tail end of the link then rotates 
into place. 
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Figure 6. Successive transformations between two links made by rotating a link so that there is a 
progressively larger angle between the links as vectors (or smaller angle made between them as lines). 
The two boundary conditions (the initial and final conditions) are shown as black links, and an intermediate 
state is shown as a red link or links. The arcs traced out by the end points are shown in blue or green, while 
straight line motions when they are not along the links themselves are shown in grey. The distance travelled 
over the course of the transformation is given below each figure. 
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In figure |7p, straight lines from the end points will not satisfy the obtuse condition in section 13.21 
so the transformation must involve rotations. Here a straight line transformation takes the link almost to 
the final state. It then must undergo a small rotation to complete the transformation. Seen in reverse, the 
vertical arrow must rotate to a critical angle determined by the criterion in section 13 .21 before the link 
can finish the transformation by pure translation. 

Figure [TJi is figure (6f once again. The final condition (the tilted link) will be systematically changed 
by translating it vertically away from the horizontal link (which we choose arbitrarily as the initial 
configuration). 

In figure [Tp the tilted link is translated a distance 1/3 vertically. The transformation can be achieved 
by rotating the horizontal link about a point equidistant from both arrowheads, to the red intermediate 
configuration. The link then rotates about the arrowhead into the final configuration. The distance is 
still the angle rotated for the reasons mentioned above in the context of figures [Sf-G, = (150/180)7r, 
which is unchanged from[7E. In fact, so long as the arrowhead can be reached by rotation (the translated 
distance is less than d where d is the solution to d 2 + d + 1 — y/3 = for this angle), then the distance 
will be unchanged. The transformation at the critical distance is shown in figure [Tfj. The rotations now 
occur about the end-points: the tail and head of the link. 

In figure |7H the translated distance is now equal to 1. The transformation first consists of a rotation 
about the tail to a critical angle (blue arc and red intermediate state), then a translation much like that in 
figure H] (grey straight lines between red intermediate states), and finally a rotation about the head (green 
arc) to the final configuration. 



4. 2-link chains 

We now consider the next simplest case of 2 links (3 beads). The Lagrangian now reads: 

£(r 1 ,r 2 ,r 3 ,r 1 ,r 2 ,r3) = ^+^+y / ^-^i2((r2-r 1 ) 2 -l)-^ 23 ((r 3 -r 2 ) 2 -l) (38) 

which has EL equations (c.f. eq.s ll5alfl"5cl) * : 

v A + l AB r B/A = (39a) 
v B - 1 AB r B/A + 1 BC r c/B = (39b) 
v c -l BC r c/B = 0. (39c) 
The corner conditions (|18al) , (fT9l) imply 

M'")=M' + ) 

so the direction of motion cannot suddenly change, unless along one part of the extremal the velocity of 
point i is zero (the point is at rest), where its direction v is then undefined. 

The boundary conditions described in section 12.3.11 hold as well, so the end points can either be at 
rest, move in straight lines, or purely rotate. This gives 3x3 = 9 possible scenarios to investigate here, 
many of which can readily be ruled out. For example consider the states in figure [8^. Because A and A' 
are in the same position, rotation and translation of A are ruled out and point A remains at rest, leaving 3 
scenarios for the other end point C. However since C and C are at different positions and ABC are along 
a straight line, C cannot remain at rest initially, leaving either translation or rotation for point C. 

Suppose C translates towards C as in figure [8J5. Then v c = and from (|39c|39b| ) 1 BC = and 

v B = l AB r B/A- B cannot move in a straight line without moving point A, so 1 AB 7^ and thus B must 
rotate about point A. The transformation then proceeds as in figure [8J5 until B reaches B' and C reaches 
C" . Then however if C" were to rotate to C, the trajectory would violate corner conditions at point C". 

X * The links have length 1 in our dimensionless formulation, so the vectors r B/A and r c/B could also have been written as unit 
vectors f B/A and f cm ■ 



Minimal distance transformations 



19 



D = 1.571 




D = 2.616 




D = 1.571 




D = 2.618 




D = 1.730 




I) = 2.616 
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Figure 7. Successive transformations between two links made by translating one link with respect to the 
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Figure 8. (a) Initial and final states for a chain of two links. The transformation in (b) is non-extremal 
because it violates a corner condition at C", (c) and (d) are degenerate minima- rotations occurring about 
B' or B both have the same length. Intermediate states is shown in red have opposite convexity in (c) and 
(d). 

Therefore the direction of translation of C must not be directly to C but must be tangential to the arc 
C'C" as in figure [8b. 

The reverse of this transformation is allowable as well, as can be seen by swapping the labels 
ABC —> A'B'C . Here C first rotates to the critical angle shown in fig [8tl and then translates to C 

In fact one can see that links BC and B'C along with lines BB' and CO form a quadrilateral as in 
figure [51 with the same consequences for rotation to a critical angle. For the links in fig [8] the situation 
is symmetric so rotation can occur at the beginning or end of the transformation. Figure [9^ shows an 
example with this symmetry broken, so that the distance is different depending where the rotation occurs, 
as in figures [5^,b. In this case, the transformation in fig [9b has the minimal distance, and that in fig [9}) is 
subminimal. Extensions of the transformation in figure [9] to large numbers of links were explored in [Q]|. 

4.1. Transformations involving a change in convexity 

Transformations between configurations with opposite convexity involve motion out of the plane, even 
if the initial and final states lie in the plane. If the transformation is constrained to lie in plane, the 
trajectories of some points will be non-monotonic- those points must move farther away from their final 
positions before approaching them. We illustrate these ideas with some examples below. 

Consider the initial and final states in figure [TOl We again imagine B rotating to B' . If C were 
to translate to C one would have the intermediate configuration A'B"C . Now C and A' must remain 
at rest to satisfy corner conditions. Then the only way to finish the transformation is for B" to rotate 
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D = 3.114 



D = 2.985 
C 



Figure 9. (a) Initial and final states for a polymer of 2 links. The angle between AB and A'B' is tt/4. The 
minimal transformations in (b) and (c) are now no longer degenerate, (c) is the global minimum. 



about the axis A'C, however then the trajectory of B violates corner conditions and so is not extremal. 



In Appendix C we take up the issue of minimal transformations for this case when the links are 
constrained to lie in a plane. 

> ^ 

We thus seek a point B" and resulting trajectory BB"B' such that arc BB" satisfies corner conditions 
with arc B"B'. 

One solution is to effectively place B" at position B' by considering the boundary condition with C at 
rest (and A at rest). Then B rotates to B' about axis AC, and the trajectory of B lies on a circle defined by 
the intercept of two unit spheres centered at A and C. The sphere about A is drawn in figure[TT]as a visual 

aid. Along arc BB' both 1 AB ^ and 1 BC ^ 0. Once in configuration A'B'C, C can then undergo rotation 
about B' to C, with A' and B' stationary. 

The transformation in [TTk is a local minimum in distance, however it is not the global minimum. 
A shorter distance transformation can be seen by considering the reverse transformation. Imagine A' 
and C stationary while B' rotates about axis A'C in figure [TTb . This rotation of B' follows a circular 
trajectory defined by the intercept of two unit spheres centered at A' and C. The rotation occurs until 



V c 




Figure 10. A transformation between two states of opposite convexity: ABC has convexity down and right, 
while A'B'C' has convexity up and left. There is no extremal transformation in the plane that can connect 
them, without some apparent violation of corner conditions. 
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point B", which is the point where above circle is tangent to a great circle on the unit sphere about A and 

passing through B. The arc BB" is a great circle because this is a geodesic for point B given A is fixed, 
which follows from the Euler equations (|39b|39c| ) when 1 BC = 0. The great circle is defined by the plane 
containing the points A, B, and B". 

The angle between the (variable) vector BC of link BC and the tangent the the arc B'B" is always tt/2, 
so once the corner condition is met, point C on link BC can move in straight line motion from C to C 
while B moves on the great circle from B" to B. That is, the quadrilaterial criterion of section I3TT1 is met 
for uBB"C'C. 

To find point B" , let its position be r B - = (x Ql y(x ),z(x )). The great circle is defined by the 

plane passing through the points A, B, and B" . This plane has normal n = AB x AB" = (1,0,0) x 
(x ,y(x ),z{x )) = (0, — z{x ),y(x )). At the point B" the normal is orthogonal to the tangent vector 
of the circle defined by rotation about the AC axis. This tangent vector is t = dr/ds = x s (l, y xi z x ) by the 
chain rule. At B", t • n = 0, or 

-z{x )y x (x ) +y{x )z x (x ) = (40) 
The functions y(x) and z(x) are defined by the intercept of two unit spheres centered at (0,0,0) and 
(1/^,1 + 1/^0), giving 



Together (gO]) and gB give 



y{*) 



V2 
2 + V2 



x 



■y{x) 



(41) 



/ y/2-l \ 
2(V2-1) 



\^2(5V2-7)J 



The distance travelled along arc BB" is Bjr , where cos BB - = x„ = \[2— 1. The distance travelled along arc 

B"B' can similarly be shown to be r B .. B - = sin(7r/8)cos _1 (2\/2 — 3). Adding the distance CC, the total 
(minimal) distance is thus V = 2.576. There is of course a degenerate solution to the above with z — > — Z. 



4.2. Transformations with initial and final states in 3-D 

We now give a representative example where the initial and final configurations do not lie in the same 
plane, as shown in figure [121 Because AB _L AA' and BC _L CC, neither A nor C will rotate about B 
as part of the transformation. Nor can ABC simultaneously translate directly to A'B'C, because for 
example quadrilateral uAA'B'B does not satisfy the rule of opposite angles > n/2, so link AB cannot 
slide (translate) to A'B'. 

This leaves 3 options for the initial stages of the transformation: 

1. ) A translates, B rotates, C remains fixed. B then rotates about C in the CBB' plane. The initial direction 
of motion of B is then v B = (— i + k)/ y2, however then v A can only move backward to preserve link 
length (v A = — k), similar to figure [ATI This rules out case (1). 

2. ) A remains fixed, B rotates, C remains fixed. B then rotates towards B' about axis AC until it reaches 
a critical angle where line B"B' is tangent to its circular trajectory (see fig. [T2k). At this point the 
quadrilateral uB"CC'B' does not have opposite obtuse angles, so a straight line transformation to A'B'C 
is not possible. It is possible to transform to a configuration A'B'C" , where C" is at position (1,1,1) and 
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Figure 11. Subminimal (a) and minimal (b) transformations for the boundary conditions in figure[T0] The 
distances for each transformation are approximately 3.007L 2 and 2.576L 2 respectively. Transformation (a) 
proceeds from ABC by first rotating B to B' about axis AC, then rotating C about point B' . Transformation 
(b) proceeds from ABC by simultaneously translating C to C' while rotating B about A on a great circle to 
point B" . Finally point B rotates from B" to B 1 about axis A'C '. 

angle ZB'C'C = ir/2, so that v c = k. Then the transformation is completed by a n/2 rotation of C" about 
B' . This transformation is subminimal. 

3.) A remains fixed, B rotates, C translates. In this case, B rotates toward B' in the BAB' plane, 
while C translates to C, until the state AB"C" is reached (see fig. \12b). State AB"C" can be found as 

follows. Because the rotation of B is about the axis (0, -1/v^, 1 /y/2), the position AB" of B" after 
rotation of the (critical) angle is (cos , sin / y/2, sin / y/2). This angle is then determined by the condition 

A5" • IFb' = 0, where b/b' = AB' -AB". The solution to this condition is simply = vr/4. The location 

> 

of C" is then determined from the condition that the link length from B" to C" is one: \B"C"\ = 1, where 

B"C" = A&' + tCC'. Solving this condition for t gives the position of C" as (^^, 1, 2(2 ~^ } ). At this 
point the quadrilateral uB'B"C"C has opposite obtuse angles, and quadrilateral uAB"B'A' has opposite 
angles = %/2, so it is in a bowtie configuration as in the end point configurations in figure|4] Therefore all 
points AB"C" can translate from this intermediate state to their final positions A'B'C. The total distance 
travelled is + \AA'\ + \CC'\ + \B"B'\ or V = 2 + vr/4 + y/5 « 5.022. The reverse of this transformation is 

also possible, where point B' rotates about A' in the plane B'AB, while C translates along C'C. Inspection 
reveals the distance covered is the same as the forward transformation. 

5. Limit of large link number 

From the transformation discussed in section |4~TI we see that if both ZABC and ZA'B'C were n/2 
as in figure [T3b . then the transformations in figures [TTk and [TTb become degenerate, having distance 
V = 7r/ y/2. The transformation is completed by a single rotation about axis 13. 

We can now examine the effect of increasing the link number. Let the number of links increase to 4, 
and let us preserve the symmetry that is present about the horizontal axis in fig [T3k . so the initial and 
final states become an octagon (figure [13b). In the limit 7Y — > oo, the figure becomes a circle. 

If we separated the links in figure [T3k by some distance in the y direction (perpendicular to axis 13), 
then the minimal transformation involves the same rotation of 2 about axis 13 until a critical angle 8 C , 
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25., 2.5-, 




Figure 12. (a) Subminimal transformation and (b) minimal transformations between ABC and A'B'C' (see 
text) 



after which all three points 123 can translate in straight lines to l'2'3'. In the same fashion, the minimal 
transformation for the octagonal transformation in fig [T3b involves a rotation of point 3 out of the plane 
about axis 24 to a critical angle 9 C at which the point is located at position 3". Once this critical angle is 
reached, point 3 translates in a straight line from 3" to 3'. 

Because points 1 and 5 are stationary to satisfy corner conditions, points 2 and 4 must move in great 
circles about points 1 and 5. However points 2 and 4 cannot finish the transformation by moving on great 
circles. At the configuration 1'2"3'4"5' in figure [T3b. point 3 has finished the transformation, but points 
2 and 4 have not. To satisfy corner conditions at the points 2" and A", the great circles must be out of 
plane as well. At points 2" and 4", the transformation finishes with rotations about axes 1'3' and 3'5'. 
The total distance V « 7.93. 

Of course the time reverse of this transformation (equivalent to swapping primed and unprimed labels) 
is also a minimal transformation, as is the transformation obtained by reflection about the z = plane. 

Now consider increasing the chain to 6 links, so the combination of r,(0) and r,(r) becomes a 
dodecagon (12-sided polygon, see figures [T3~b-d). As before the midpoint vertex (here r 4 ) must rotate out 
of the plane about axis 35 to a critical angle c before translating in a straight line to r 4 /. This critical angle 
is where 34" • 4 // 4 / = 54" • 4"4' = 0. The quadrilaterals □22'3'3 and c^SS^' are of the type in figure [5J 
so point 3 must rotate about r 2 (0) to a critical angle where 23" • 3"3' = 0, and likewise for point 5. 

While point 3 rotates to its critical angle, point 4 translates along line 4"4'. Points ri(0) and r 7 (0) 
overlap with ri (T) and r 7 (r) and so remain fixed to satisfy corner conditions. After point 3 has reached 
its critical angle, it can translate along 3"3' as point 2 rotates about ri. However to satisfy corner 
conditions at point 2", the rotation cannot remain in the x — y plane. Point r 2 » is determined as the 

point where t • n p i ane = 0, where t is the tangent to the arc 22" defined by rotation about axis 13', and 
Wpiane is the normal to the plane 122", i.e. x r 2 »/i. The same process holds for point 6. These critical 
points and some intermediate states for the transformation are shown in figure [T3ti . The total distance 
covered by the transformation is V m 16.3. 

It is sensible to consider the total length of chain as fixed to say L = 1, and to let the link length ds N 
for the chain of /V links be determined by Nds N = L. Because distances scale as ds^, the /V = 2, 4, 6 cases 
have T>2 ~ 0.555L 2 , V4 « 0.496L 2 , m 0.445L 2 . Note that this distance decreases with increasing 
number of links: the constraints on the motion of the various beads during the transformation are relaxed 
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Figure 13. Examples of transformations between initial and final states of opposite convexity, for increasing 
numbers of links, (a) illustrates the transformation for N = 2 links. (b)N = 4 and initial and final state form 
an octagon. (c,d) N = 6 and initial and final states form a dodecagon, (c) top view, (d) view in perspective. 
Rotations are shown as solid color lines (either green or blue). Translations are shown as dashed lines. The 
grey dashed lines underneath 3"3' in (b) and 4"4' in (d) are shown only to illustrate that those lines are 
above the plane. 

as the number of links is increased. 

We can then imagine resting a piece of string on a table in the shape of a semi-circular arc, and then 
asking how one can move this string to a facing semicircle of opposite convexity. So long as the string 
has some non-zero persistence length £ P , the transformation of minimal distance must involve lifting the 
string off of the table to change its local convexity. The vertical height the string must be lifted (see 
fig[T3tr) is of order ~ sm(%£ P /L) ~ £ P /L, which goes to zero for an infinitely long chain. 

As the number of links N — > oo, some simplifications emerge. In particular the contribution to the 
total distance due to rotations becomes negligible, and the translational component dominates. To see 
this note that the distance due to straight line motion scales as: 

£>(st. line) ~dsNL~L 2 

while the distance travelled during rotations scales as 

X?(rot.) ~ dsN(-ds) ~ L 2 /N 

where we assume the worst case scenario where an extensive number of links must rotate before 
translating. Because translation dominates the distance as N — > oo, the distance travelled converges 
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to L times the mean root square distance (MRSD), i.e. 



N+l 



i=l 



= L ^Ev /(rB '- rA ' )2 



L(MRSD) 



(42) 



The MRSD for the examples in figures [T3b.d are 0.394L and 0.400L respectively, which are both 
less than the actual distances travelled (in units of L). In the limit N — > oo, where the polygon becomes 
a circle, the distance converges to = 4L 2 /% 2 »s 0.4053L 2 . For large N systems then, it is a good first 
approximation to use MRSD for the distance. 

The MRSD is always less than the root mean square distance (RMSD), except in special cases when 
they are equal. To see this, we can apply Holder's inequality 



where gk,hk ^ 0, a, (3 > 0, and a -\- (3 — 1. With the specific identifications gk — (r B £ — v A k ) 2 = Art, 
h k =l, and a = (3 = 1/2, we have directly 



For example the RMSD for the circle configuration discussed above is \flLjii « 0.4502L, which is 
greater than the MRSD. 

The fact that the distance converges for large N to MRSD rather than RMSD suggests that RMSD may 
not be the best metric for determining similarity between molecular structures, although it is ubiquitously 
used. This fact warrants future investigation- it has implications in research areas from structural 
alignment based pharmacophore identification (3l|4l|5|| to protein structure and function prediction 

It was shown in flU that chains with persistence length characterized by some radius of curvature 
R have extensive corrections to the MRSD-derived minimal distance, which do not vanish as N — > oo, 
but remain so long as R/L is nonzero. Likewise, chains that cannot cross themselves have nonlocal EL 
equations and extensive corrections to the minimal distance. Nevertheless, it is worthwhile to investigate 
some more complex polymers with MRSD as an approximate distance metric. We pursue this in the next 
section. 

5.1. MRSD as a metric for protein folding 

Here we examine the use of MRSD as a metric or order parameter for protein folding. To this end 
we adopt an unfrustrated C a model of segment 84 — 140 of src tyro sine-protein kinease (src-SH3), by 
applying a Go-like Hamiltonian [|8l |9l Q23 to an off-lattice coarse-grained representation of the src-SH3 
native structure (pdb lfmk). Amino acids are represented as single beads centered at their C a positions. 
The Go-like energy of a protein configuration a is given by the following Hamiltonian, which we will 
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explain term by term: 

H(a\N) = 

+ 



+ 

Adjacent beads are strung together into a polymer through harmonic bond interactions that preserve 
native bond distances between consecutive C a residues. Here r a and r N represent the distances between 
two subsequent residues in configurations a and the native state N. As with other parameters in the 
Hamiltonian, the distances r N are based on the pdb structure and may vary pair to pair. The angles N 
represent the angles formed by three subsequent C a residues in the pdb structure, and the angles cj) N 
represent the dihedral angles defined by four subsequent residues. The dihedral potential consists of a 
sum of two terms, one with period 2ix and another with 2tt/3, which give cis and trans conformations 
for angles between successive planes of three amino acids, with a global dihedral potential minimum at 

(p N E [-7T,7T]. 

The parameters k r , kg, and k^, are taken to accurately describe the energetics of the protein backbone: 

we used the values k r = 50kcal/mol, kg = 20kcal/mol, k^ = 1 kcal/mol and k^ = 0.5kcal/mol for 
molecular dynamics (MD) simulations using the AMBER software package. For MD simulations using 
LAMMPS, we had used slightly different values: k r = 80 kcal/mol, kg = 16kcal/mol, k^' = 0.8kcal/mol 

and jfc® =0.4 kcal/mol. 

<p 

The last line in equation (1431 ) deals with non-local interactions, both native and non-native. If two 
amino acids are separated by 3 more along the chain (\i — j\ > 3), and have one or more pairs of heavy 
atoms within a cut-off distance of r c = 4.8 A in the pdb structure, the amino acids are said to have a native 
contact. Then the respective coarse-grained C a residues are given a Lennard-Jones-like 10-12 potential 
of depth e N = —0.6 kcal/mol (—0.8 kcal/mol for LAMMPS simulations) and a position of the potential 
minimum equal to the distance of the C a atoms in the pdb structure. That is, cr, 7 is taken equal to native 
distance between C a residues i and j if i- j have a native contact. 

If two amino acids are not in contact, their respective C a residues sterically repel each other 
(cnn = +0.6 kcal/mol). Thus e NN = if i- j is a native residue pair, while e N = if i- j is a non-native 
pair. For non-native residue pairs, ay = 4 Angstroms. 

In an arbitrary configuration a, two C a residues i and j are considered to have formed a native contact 
if they have a distance ry < 1.2ay. The results do not strongly depend on the specific value of this cutoff. 
The fraction of native contacts present in the particular configuration a is then defined as Q (or Q a ). 

The MRSD of configuration a is found by aligning this configuration to the native structure, by 
minimizing MRSD over 3 translational and 3 rotational degrees of freedom. 

Constant temperature molecular dynamics simulations were run for this system using both AMBER 
and LAMMPS simulation packages. The probability for the system to have given values of Q and 
MRSD within (£>,£> + AQ) and (MRSD, MRSD + AMRSD) is proportional to the exponential of the free 
energy F(Q,MRSD). Thus the free energy can be directly obtained by sampling, binning, and taking the 
logarithm: 

F(Qi,MRSDi) -F(Q 2 ,MRSD 2 ) = -k B Tlog ( P ^Zlln\ ) (44) 

\p{Q 2l MRSD2)J 

with F(1,0) = E N , the energy of the native structure. 

Figure [14] shows the free energy surfaces obtained using the above recipe, for the AMBER (fig [Ha) 
and LAMMPS (figfT4b) molecular dynamics routines. The temperature is taken to be the transition or 
folding temperature T F , where the unfolded and folded free energies are equal. 
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Notice that F(Q) is comparable for both as it should be, moreover F(MRSD) is as well. However the 
free energy surface plotted as a function of both Q and MRSD shows a marked difference. In addition 
to a native minimum, the LAMMPS routine has an additional minimum at Q « 0.95 and MRSD « 8.4. 
The conformational states in this bin are closely related, with an average MRSD between them of 1.8A. 
We can take the most representative state in this bin as that which has a minimum MRSD from all the 

others in the bin (at Q « .95, MRSD « 8.4): min (^'^MRSDij/ Y?#i) ~ L6A - Inspection reveals that 
this state is a mirror image of the pdb structure (see fig fT4b): If we reflect this structure about one plane, 
and subsequently align this reflected structure to the pdb one, the MRSD is only 1.1 A. 

The discrepancy in free energy surfaces corresponding to the presence of a low energy mirror-image 
structure arises because the COMPASS class 2 dihedral potentials in the LAMMPS algorithm do not 
ascribe a sign to the angle (ft, so the full range [— 7r,7r], is projected onto [0,7r]. This gives the set of 
actual dihedral angles {</>,• + 7r} the same energy as the set {(pi}, so that the dihedral potentials have two 
minima rather than one, and thus a protein chain of the opposite chirality (a mirror image) is allowed and 
has the same energy as the pdb structure. We found that the CHARMM and harmonic dihedral styles do 
not have this problem, however they have less versatile function forms, so that we favored modifying the 
COMPASS dihedrals to define <fi over its full range. 

6. Conclusions 

Analogously to the distance between two points, the distance between two finite length space curves is a 
variational problem, and may be calculated by minimizing a functional of 2 independent variables s and 
t, where s is the arc-length along the chain, and t is the 'elapsed time' during the transformation. 

We derived the Euler-Lagrange (EL) equation giving the solution to this problem, which is a vector 
partial differential equation, with extremal solution r*(s,t). We also derived the sufficient conditions for 
the extremal solution to be a minimum, through the Jacobi equation. Once the minimal transformation 
r*(s,t) is known, the distance V* = V[r*] follows. 

We provided a general recipe for the solution to the EL equation using the method of lines. The 
resulting N+l EL equations for the discretized chain are ODEs that can be interpreted geometrically 
and solved for minimal solutions. Solutions consist generally of rotations and translations pieced together 
so the direction of velocity of any link end point does not suddenly change (the Weierstrass-Erdmann 
corner conditions). 

We explored the minimal transformations for the simplest polymers, consisting of 1 or 2 links, in 
depth. For transformations between 2 links, convexity becomes an issue (the analog to the direction of 
the radius of curvature for a continuous string). For example, even if the initial and final states lie in the 
same plane, if the convexities of these states are of opposite sign the transformation must pass through 
intermediate states that are out of the plane. Similarly, given a semicircular piece of string lying on a 
table, to move it to a semicircle of opposite convexity using the minimal amount of motion, the string 
must be lifted off the table. 

The study of minimal transformations between small numbers of links has applications to the inverse 
kinematic problem in robotics and movement control. In the inverse kinematic problem, one is given 
the initial and final positions of the end-effector (the hand of the robot), and asked for the functional 
form of the joint variables for all intermediate states. Generally there is no unique solution until some 
optimization functional is introduced, such as minimizing the time rate of change of acceleration (the 
jerk), torque, or muscle tension (see the review [11J and references therein). The minimal distance 
transformation would be relevant if one sought the fastest transformation between initial and final states, 
without explicit regard to mechanical limitations. The indeterminate intermediate points can be handled 
variationally as a free boundary value problem. 

In the limit of a large number of links, some simplifications emerge. For chains without curvature 
or non-crossing constraints, the distance converges to L times the mean root square distance (MRSD) of 
the initial and final conformations. So for example the distance between 2 strings of length L forming 
the top and bottom halves of a circle respectively is 4L 2 /n 2 , the distance between horizontal and vertical 
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Figure 14. Free energy surfaces for the folding of Go-model src-SH3 using two molecular dynamics 
simulation packages, AMBER (a) and LAMMPS (b). The contour plots give F(Q,MRSD). The projections 
F(Q) and F(MRSD) are also shown on each side. The COMPASS class 2 dihedral potential in LAMMPS 
allows for a mirror image of the folded structure (red color structure in inset) that is not immediately evident 
from the F(Q) or F(MRSD) surfaces. Future implementations of LAMMPS using COMPASS dihedrals 
for biomolecular simulations must then correct for dihedral angles defined on the interval [— w, tt]. 
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straight lines of length L which touch at one end is L 2 / \fl, and the distance to fold a straight line upon 
itself (to form a hairpin) is L 2 /4. 

The fact that for large N the distance (over L) converges to MRSD rather than RMSD suggests that 
RMSD may not be the best metric for determining similarity between molecular structures, although it 
is ubiquitously used. Adopting MRSD may lead to improvements in structural alignment algorithms. 

The MRSD was investigated as an approximate metric for protein folding. Free energy surfaces 
for folding were constructed for two simulation packages, AMBER and LAMMPS. It was found that 
including MRSD as an order parameter uncovered discrepancies between the two molecular dynamics 
algorithms. Because dihedral angles in LAMMPS (at least in COMPASS class 2 style) are only defined 
on [0,7r], the potential admits a mirror image structure degenerate in energy with the native structure. 
This is easily remedied and should not be interpreted as a deficiency in the LAMMPS simulation package 
so long as one is aware of it. It should be mentioned that the mirror-image structure would also have 
been seen had RMSD been used as an additional order parameter. 

It will be important for future studies to address the effects of persistence length and non-crossing on 
the distance between biopolymer conformations JT]. Also important is the role of entropy of paths or 
transformations in describing the accessibility of a particular biomolecular structure. Along these lines it 
will be interesting to investigate whether the distance can be a predictor of folding kinetics, or proximity 
to the native structure. 

It is also an interesting question to ask whether the actual dynamics between polymer configuraitons 
resembles the minimal transformation, after a suitable averaging over trajectories. This question is linked 
with the role of the entropy of transformations described above. It is also related to the problem of finding 
the dominant pathway for a chemical reaction IfTZll . which has recently been applied to the problem 
of protein folding [fL3l . We have focused here on the question of geometrical distance for complex 
systems, which can be separated from the calculation of quantities such as reaction paths that depend 
intrinsically on energetics, i.e. on the specific Hamiltonian of the system. Quantifying the relationship 
between geometrical distance and the dominant reaction path is an interesting future question worthy of 
investigation. 

The notion of distance and corresponding optimal transformation for a system with many degrees of 
freedom is fundamental to a diverse array of research subjects. Hence we saw potential applications 
for this metric in areas ranging from drug design to robotics. It is not clear at present how useful 
the calculation of the true Euclidean distance between high-dimensional objects will be for practical 
applications, but we are optimistic. 
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Appendix A. Necessary conditions for straight line transformations 

It was shown in section 137X1 that to have straight line transformations between links, it is sufficient to have 
facing obtuse angles on opposite sides of the the quadrilateral defined by the transformation as shown in 
figure |3K- We now show that it is a necessary condition as well, i.e. we show that a slide in the correct 
direction is not possible in the absence of obtuse angles. 

Without loss of generality assume that the link is initially along the z axis. The paths travelled by the 
link ends are shown in the figure. Note that the end point trajectories of A and B are in 3D space so the 
paths travelled by A and B need not cross or lie in the same plane. Let the unit vector along A's path 
be v A and the unit vector along 5's path be v„. Because the angles that the path of A and the path of B 
make with the link are acute, the z-component of v B (= z B ) is negative and the z-component of v A (z A ) is 
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Figure Al. 

positive. One can write v A and v B as 

v A = p A + z A i 
v B = p B + z B z 

where p A and p B are vectors in xy plane and z A > and z B < 0. 

Let r A (f) and r B (f) denote the positions of the A and B ends at time t: 

r A = t\ A 
r B = g (t)\ B + z 
The rigid link constraint dictates that 

(r A -r B ) • (r A -r B ) = 1 

which translates to: 

g 2 + 2g(z B -t{c + z A z B )) -2tz A + t 2 +l = 1 
with c = p A - p B . Solving for g as a function of t, keeping in mind that g(0) = 0: 



(z B -t(c + z A z B )) +\/(z B -t(c + z A z B )) 2 -t 2 + 2tz A . 

Now if g'(t) > it means that the B-end of the link is travelling in the assumed direction, and if 
g\t) < it means that B-end is travelling in the opposite direction (which means that the angle is not 
acute anymore). Writing g'(0) we get: 

//r>\ 2z B c + 2z A z B 2z A Za ~ 

g(0) = — r + c + z A Z B = —-<0. 

2 Kb | Kb] 

Thus point B can only travel in the opposite direction from what was assumed, which in turn means 
an all-acute slide is not possible. We conclude that the condition of "facing obtuse angles" is necessary 
and sufficient for transformations consisting only of pure translations. 



Appendix B. Critical angles 

The concept of critical angle was first introduced in 13 .21 In order for a straight-line slide of both ends 
to be possible, at some stage during the transformation the link needs to rotate about one of the ends, 
with the other end being stationary. In principle the rotation can be about either of the two ends and it 
can happen at the beginning or the end of the transformation. The conditions on the critical angle or 
orientation can be readily derived from the broken extremal conditions. It was seen from ll8a| and[T9l the 
non-trivial corner conditions read: 

v/| + =v«l_. (B.l) 
We know that the path travelled by the moving bead during the rotation is circular and the path that is 
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Figure Bl. Transformation in which both ends stay on a linear track 

travelled during the slide part is a straight line. Broken extremal condition forces these two paths to be 
patched smoothly, which means that the straight-line path should be tangent to the circle. In the 3D case, 
for the broken extremal condition to be satisfied, the straight line slide path and the circular rotation path 
should lie in the same plane. For example in figure [5] where B is rotating about A initially to B\ and then 
slides to B', the rotation has to be in the plane formed by the three points ABB'. 

Matching the directions of velocity as in (|B.1I) does not itself mean that a link can subsequently 
slide in a straight line, however at the tangent point, the tangent line to the circle is perpendicular to 
the radius, hence one satisfies this second condition as well. Below we derive an analytical expression 
for the critical angle for a particular case of single link problem, as an example and illustration of the 



discussed concepts. Furthermore the particular example will be used later in Appendix C to introduce 
minimal transformations in 2 dimensions. 

Consider the single link action with the particular parametrization s = s(9) , as discussed in section 




s 2 + 1 + IscosQ + Vs 2 ) dO. (B.2) 



where s = A()A is the (signed) distance of A-end from its initial position, and 6 is the angle between the 
link and the horizontal line (see figure IBTI) . 

The Euler Lagrange equation of motion reads: 

d , s s'+ cosQ 

d9 + 7? + 1 + IscosB ~ 
We consider a transformation which is not (necessarily) a minimum: 

s = acos#-sin# + Z? (B.4) 

with a and b parameters to be determined. 

Such a transformation in fact forces the two ends to travel on a straight line (right from the beginning) , 
but the A side may in fact retreat and then move forward. We call such a transformation a "hyperextended 
transformation". A sample transformation of this kind is shown in figure IBTI The parameters a and b in 
(IB .41) can be tuned to meet the boundary conditions (see below). 

In fact it is seen that point A on the link retreats backwards until it reaches some critical angle, which 
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is when link AB makes an angle f with the straight line BB' that point B travels on. Subsequently A then 
moves forward towards A'. 

Assume that 9 runs from 0i to 2 , where < 9 2 < it/2. For simplicity assume that both these angles 
are between and ~ . 

The boundary conditions dictate that: 

s(6 l )=0 (B.5) 

s{9 2 ) = I (B.6) 

where / is the distance between A and A', 
a and b can be explicitly solved to give: 

— sin 2 + sin Q\ — I 



cos 9 1 — cos 6*2 



(B.7) 



cos 6\ (—sin 6? — I) + sin6>i cos 8 2 

b= 2i a ' B -^ 

COS fc^i — COS 02 

For our purposes we only need to note that the critical angle occurs when s = ^ becomes zero, that 
is when A stops going backward and starts moving forward: 

s = -a sin -cos = (B.9) 

where a is given in lB.7[ 

We can now ask what should 9\ be so that there is no need for the link to go backward, i.e. it moves 
forward from the beginning and the transformation is monotonic. Equations (IB .91) and (IB.7I) give: 

— sin 2 + sin0 — / 
cos0+ 1 sin0 = O (B.10) 

COS — COS 2 




Figure B2. Geometric proof for critical angle condition 
For pedagogical reasons we prove condition (IB. 101) using analytic geometry as well. Looking at 
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figure IB21 we have the following: 

g 2 + l\ = l (B.ll) 
g 2 + l 2 = a 2 (B.12) 

a I + I1+I2 



We can solve g = a/1 — 1\ and a = y/l — lf + l^ from the first two equations and substitute in the third 
equation to give: 

/= Vi-g+*2 _r h (B14) 

On the other hand based on our results for g and a we have: 

/1 _/2 

sinfli = — ^= L= (B.15) 

cosfli = — = J? (B.16) 

sin 2 =/i (B.17) 



cos0 2 = yl-^ (B.18) 

Substitution of eqns (|B.15HB.18b in equation (IB. 101) gives equation (IB. 141) after some simplification. 

For the particular case that we have discussed, the proposed transformation is in fact a minimal 
solution if 9\ is greater than the critical angle, because in that case a simple slide would be possible. If Q\ 
is less than the critical angle a locally minimum solution as we know is pure rotation to the critical angle 
and then straight line slide. Pure rotation has a nice geometric interpretation in our parametrization. it 
corresponds to the null solution s = 0. Since at the critical angle s = we see that s = will be smoothly 
patched with s = a cos 9 — sinO + b, as mandated by the corner conditions in equation (I18al) . 



Critical Angle 

Figure B3. A minimal transformation in s(0) parametrization. The horizontal segment corresponds to pure 
rotation and the curved section corresponds to slide on straight paths. Here the corner conditions demand 
that the derivative s be continuous at the critical angle. 



Appendix C. Minimal transformations in 2 dimensions 

It was seen in section 14.11 that for the case of two links when one is confined to moving in a plane, 
satisfying the constant link length constraints and corner conditions do not seem to lead to solutions 
which are extremal. However given the additional constraint that the links must lie in a plane, there must 
be one or a set of minimal transformations. We need to look at other forms of transformations, namely 
compound straight line transformations. We will elaborate on the idea starting with single links. 
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The hyper extended solution that was discussed previously in | Appendix B| can be considered as a 
very special example of compound straight line transformation. These are transformations that are made 
strictly from straight line paths with no pure rotation. A more general transformation is shown figure 
beside the old transformation. 




Figure C 1 : The previous hyper extended solution is 
shown along with a more general compound straight-line transformation, where AA" travels in some general 
direction. Length of each line segment is written beside it. For the hyper extended solution the value of AA" is 

multiplied by two because the path is travelled twice. 

Note that the corners do not technically violate the corner conditions because the speed of "A" bead 
is zero at the corner point in any parametrization that can simultaneously describe A motion and B 
motion: Since at the corner point, the link makes an angle of 90 degrees with the path that B travels, the 
speed of B at the critical angle in infinitely larger than the speed of A. In fact one sees that we have an 
instantaneous pure rotation about A-bead, when it is at the corner point. v a is not clearly defined at the 
corners, and everywhere else (when the speed of the bead(s) is not zero), the two beads are travelling on 
a straight line. The two solutions depicted in the figure come from two different parametrization of the 
most general form of the action and result in different distances. But each of them is a local minimum 

once the direction of AA" is picked, and these local minima have different values for the distance. 

We can then ask about the best position to put the corner point, to minimize the distance travelled in 
the compound straight line transformation, with respect to other compound straight line transformations. 
We assume the corner occurs on one side and we take it to be the "A" side. 

Note that at the corner, the link makes a 90° angle with the B-bead path BB', meaning that the distance 
from the corner point to B path is always the length of the link, i.e. unity here. Also note that the total 
distance that the "A" -bead travels is the distance from the initial point A to the corner point A", plus the 
distance from A" to the final position A'. 

The locus of points with equal sum of distances from two points A and A' defines an ellipse with foci 
at A and A'. Moreover the length of the major axis of the ellipse equals the sum of the distances from 
the foci. Thus the smaller the major axis of the ellipse with foci A and A', the smaller the total distance 
travelled by the "A" -bead. Moreover A" should sit on a line parallel to B-path at a distance of 1 from the 
B-path line BB'. So in seeking the shortest distance travelled the A end of the link, we seek the point A" 
such that it lies on an ellipse with foci A and A', the ellipse shares at least one point with a line parallel 
to BB' and distance 1 away from it, and lastly that the ellipse has the smallest possible major axis (see 
figure IC2l) . So the ellipse giving the minimal distance is tangent to the parallel line, and A" is the tangent 
point. This is illustrated in figure |C2| 

This solution can be straightforwardly extended to 2 links, as depicted in figure IC3l Consider then 
the example in figure [T3k, where the links are no longer allowed to move out of the plane (see figure |C4l) . 
Here r A = r A i and r c = r C / and the above ellipses turns into a circles centered at A and C. The circles 
have radii 1 — 1 / y/2, so that the perpendicular distance from line BB' to the farthest point on the circle is 
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B' 




Figure C2. Optimal Compound Straight Line Transformation 
1 and a fully extended intermediate state is allowed. 




Figure C3. An optimal compound Straight line solution for 2 link. For this particular class of solutions, 
the problem is divided into to disjoint problems (one for each link) and solved separately. 
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Figure C4. Minimal transformation restricted to 2 dimensions, for 2 links of opposite convexity which 
form opposite sides of a square. 
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